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摘 要 : 中 国 是 世界 上 受 


早 影响 最 严重 的 国家 之 一 ,干旱 频 发 给 我 国 经 


济 社会 发 展 和 生态 环境 造成 严重 影响 。 


为 分 析 近 40 a 干旱 事件 的 时 空 特征 


,本 文 结合 三 维 聚 类 算法 ,从 
2020 年 间 干旱 事件 并 定量 分 析 干 时 事件 的 时 空 动态 演变 过 程 。 主 要 结论 如 下 :三 维 聚 类 外 


量 事件 时 空 联动 的 本 质 出 发 ,识别 中 国 1981 一 
法 能 有 效 识 别 干旱 事件 


及 其 动态 变化 过 程 。 中 国 1981 一 2020 年 间 发 生 持续 2 个 月 及 以 上 的 干旱 事件 102 场 ,空间 上 ,干旱 事件 空间 轨迹 倾 
向 于 自 东 向 西 发 展 ;时间 上 ,干旱 事件 时 间 重 蔷 度 较 高 ,长 历时 干旱 多 具有 多 峰 特 点 。 此 外 ,覆盖 范围 广 旦 严重 度 
高 的 干旱 事件 集中 发 生 于 2005 一 2010 年 。 本 文 结论 有 助 于 发 现 中 国 干旱 事件 的 时 空 演 化 规律 ,为 我 国 干旱 监测 和 


干旱 风险 管理 提供 科学 参考 。 
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政府 间 气 候 变化 专门 委员 会 (Intergovernmental 
Panel on Climate Change, IPCC ) 将 干旱 定义 为 “一段 
时 间 是 以 导致 严重 水 文 失 衡 的 异常 干旱 天 气 ”"。 
近年 来 , 随 着 全 球 气温 的 升 高 ,干旱 事件 发 生 的 频 
次 ,覆盖 面积 均 有 增加 ”。 中 国 地 人 处 典型 季风 气候 
区 ,干旱 灾害 的 影响 突出 ,每 年 造成 的 粮食 减产 从 
数 百 万 吨 到 数 千 万 吨 不 止 ,干旱 造成 的 直接 经 济 损 
失 高 达 440x10 元 'a，…” 。 科 学 识别 干旱 事件 并 定 
量 描述 其 特征 是 分 析 干 旱 时 空 演变 规律 的 前 提 , 能 
否 准确 识别 干旱 事件 ,并 定量 化 分 析 干 旱 事 件 的 时 
空 特征 及 其 发 展 过 程 是 干旱 防 灾 减 灾 和 预测 的 前 
ЖЕ 

В Rs BU CROKOCHUR , — ATA 
常 被 称 为 一 次 干旱 事件 。 要 深化 干旱 变化 规律 的 
认识 ,首先 需 准确 识别 并 提取 干旱 事件 ,而 后 进行 
多 特征 分 析 。 一 场 干 旱 事件 实质 上 是 一 个 时 空 联 
动 的 过 程 , 既 在 时 间 维 度 上 演化 也 在 空间 维度 上 变 
化 。 传 统 干旱 事件 识别 方法 多 样 ,比如 闽 值 水 平 
法 ,游程 理论 法 ,小波 变 换 法 经验 正 交 隐 数 等 。 大 
量 学 者 基于 以 上 方法 开展 了 干旱 监测 的 相关 研究 ， 
为 干旱 事件 识别 和 干旱 特征 定量 化 分 析 提 供 了 重 
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要 参考 ,如 任 天 霞 等 “基于 百 分 位 国 值 法 分 析 西 北 
干旱 区 极端 气候 变化 趋势 ; 李 天 水 等 "基于 游程 理 
论 和 Copula 函数 的 联合 应 用 ,为 区 域 干 旱 风险 评估 
提供 参考 ; Huang 等 利用 交叉 小 波 分 析 方 法 对 渭 
河流 域 水 文王 旱 与 气象 干旱 的 相关 性 进行 了 研 
究 。 经 验 正 交 函 数 (EOF) 和 旋转 经 验 正 交 冰 数 (RE- 
OF) 也 被 广泛 应 用 于 气象 要 素 时 空 特征 分 析 ””。 
然而 ,以 上 方法 局 限于 一 维 或 二 维 层 面 ,多 通过 将 
三 维 空间 压缩 为 二 维 或 一 维 来 简化 干旱 事件 的 三 
维 结构 。 比 如 固定 区 域 干 旱 指数 的 时 间 变 化 (一 
维 ) 或 某 一 时 间 有 段 内 干 湿 状况 的 空间 分 布 图 (二 维 ) 
等 。 而 干旱 事件 是 一 个 三 维 尺 度 下 的 干 湿 异 常 问 
题 , 既 在 空间 上 发 展 也 在 时 间 上 变化 , 降 维 处 理 固 
然 能 够 更 为 迅速 地 提取 干旱 事件 并 进行 分 析 , 但 同 
时 也 会 破坏 干旱 事件 的 三 维 时 空 结构 ,进而 导致 干 
旱 路 径 .发展 过 程 移动 方向 等 关键 信息 的 丢失 。 
为 解决 上 述 问题 ,近年 来 ,部 分 国内 外 学 者 聚 
焦 三 维 视角 下 的 干旱 事件 识别 与 特征 定量 化 研究 ， 
并 且 取 得 了 一 定 的 进展 。Andreadis 59" 8| HR% 
算法 识别 干旱 事件 及 其 空间 范围 ,并 建立 了 干旱 严 
重 程度 -面积 -持续 时 间 (SAD) 曲 线 ,将 每 一 次 干旱 
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的 面积 与 其 严重 程度 联系 起 来 。 后 来 ,Lloyd- 


近 减 少 到 800 mm, 大兴安岭- 阴山- 兰州- 青藏 高 原 


Hughes 等 "扩展 了 SAD 干旱 识别 方法 ,明确 了 干旱 
事件 的 三 维 (经 度 .纬度 时间) 结构 ,并 通过 逐一 识 
别 干 旱 指数 三 维 立方 体 (3x3x3) 内 的 干旱 像 元 点 的 
方式 实现 干旱 事件 的 识别 ;Cuo 等 “对 传统 三 维 聚 
类 算法 进行 了 改进 ,从 历时 严重 度 .烈度 .质心 等 


东南 部 地 区 降水 量 约 为 400 mm, 贺兰山 -青藏 高 原 
中 部 200 mm ,西海 降水 量 减少 到 200 mm ,而 塔里木 
盆地 年 降水 量 则 不 足 50 mm。 早 小 灾害 呈 渐 进 变 
化 ,分 布 在 全 国 各 地 ,灾害 从 北向 南 更 加 多 样 化 。 
干旱 是 中 国 北方 地 区 的 主要 灾害 类 型 ,而 洪水 是 南 


多 个 变量 来 度量 干旱 事件 ,并 分 析 了 亚洲 中 部 地 区 
干旱 事件 的 时 空 演 化 特征 ;Liu 等 "通过 三 维 空间 
内 气象 干旱 聚 类 和 水 文 干旱 聚 类 的 重 受 情况 研究 
了 气象 干旱 向 水 文 干旱 的 演变 过 程 ;Wen 等 “绘制 
三 维 阔 值 等 值 面 , 聚 类 时 空 互联 网 格 点 ,形成 三 维 
干旱 结构 ,表征 干旱 事件 ,描述 了 干旱 发 展 过 程 中 
的 时 空 分 布 格局 和 时 空 变异 性 ;Xu 等 "基于 三 维 聚 
类 方法 识别 了 1961—2012 年 中 国 季风 区 的 干旱 事 
件 ,并 通过 历时 、 受 影响 面积 .严重 程度 .强度 和 质 
心 等 多 个 特征 变量 来 度量 干旱 事件 ; 冯 凯 等 "基于 
三 维和 干旱 识别 方法 分 析 了 黑河 流域 干旱 时 空 连续 
动态 演变 特征 ; 邓 浴 玲 等 ”将 图 像 三 维 连通 性 识别 
方法 应 用 于 长 江 流域 的 气象 干旱 事件 的 识别 和 
提取 。 

随 着 气候 变化 和 人 类 活动 的 不 断 加 剧 ,中 国 的 
于 旱 事件 不 断 频 发 且 有 进一步 加 剧 的 趋势 ,干旱 时 
空 特征 日 趋 复杂 ,然而 中 国 干旱 事件 时 空 演 变 规 律 
仍 不 清晰 ,对 干旱 发 生发 展 过 程 的 认识 仍 有 待 深 
入 。 因 此 ,本 文 基于 三 维 聚 类 算法 从 时 空 联动 的 视 
角 出 发 识别 中 国 1981 一 2020 年 间 的 干旱 事件 ,在 分 
析 干 星 发 生 时 间 、 历 时 、 严 重度 烈度 等 常规 静态 特 
征 的 基础 上 ,定量 化 研究 干旱 事件 的 发 展 过 程 移 
动 路 径 等 动态 演变 特征 , 旨 在 深化 对 干旱 事件 时 空 
演变 规律 的 认识 ,促进 干旱 识别 与 特征 定量 化 的 发 
展 , 为 干旱 风险 管理 和 和 干旱 预测 提供 科学 参考 。 


1 研究 区 与 数据 


1.1 研究 区 概况 

中 国 位 于 亚洲 东部 .太平 详 西岸 ,地 理 坐 标 介 
F 73933'-135?05'E,3951'-53?33'N 之 间 ,研究 区 域 
不 包括 台湾 。 地 势 西高 东 低 , 呈 阶 梯 状 分 布 ,高 程 
空间 差异 大 。 中 国 主要 以 季风 气候 为 主 , 南 北 跨 纬 
度 广 ,太阳 辐射 能 量 空间 异 质 性 强 。 由 于 海陆 热力 
性 质 差 异 ,降水 量 从 东南 到 西北 逐渐 降低 ,东南 沿 
海地 区 年 降水 量 在 1600 mm 以 上 ,在 秦岭 -淮河 附 


方 地 区 最 重要 的 灾害 类 型 "5 。 半 个 世纪 以 来 ,我 国 
高 温 、 强 降水 等 极端 天 气 气候 事件 趋 多 、 趋 强 "”, 即 
使 在 一 些 通常 被 认为 潮湿 的 地 区 干旱 灾害 现象 也 
MERE. 干旱 事件 通常 集中 发 生 在 华北 平原 
至 长 江 下 游 地 区 ,华北 平原 西部 .黄土 高 原 . 四 川 盆 
地 和 云贵 高 原 有 明显 的 干旱 趋势 汪 。 频 发 的 干旱 
带 来 了 严重 的 损失 ,例如 2022 年 夏季 长 江 流域 发 生 
的 罕见 的 重大 干旱 事件 ,对 当地 工农 业 生 产 、 居 民 
生活 .生态 安全 等 造成 严重 影响 ”2 。 

1.2 数据 来 源 与 处 理 

1.2.1 中 国 逐日 降水 分 析 产 品 CPAP ФЕН 
水 分 析 产 品 (CPAP) 来 源 于 中 国 气象 数据 网 (www-. 
data.cma.cn) , 2 [B] 43913 0.25°х0.25°. CPAP 由 国 
家 气象 信息 中 心 (NMIC) 和 中 国 气 象 局 (CMA ) 基 于 
约 2400 个 气象 站 点 的 降水 数据 ,结合 最 优 插值 技术 
(О) 生产 得 到 。 图 1 给 出 了 中 国 地 区 用 于 生成 
CPAP 的 气象 站 点 分 布 。CPAP 所 用 气象 资料 均 经 
过 极端 值 的 检验 ,内 部 一 致 性 检验 和 空间 一 致 性 检 
验 的 三 级 质量 控制 吕 。 经 验证 与 独立 规范 观测 相 
FE , CPAP 1E 0.59 BE F ÉSEDGE 2S (RB) y 3.21967" 
WKH , CPAP < Н ӨК ЖЕЗ P Bsp [e i ІХ fn] 
观测 结果 有 很 好 的 一 致 性 所 。 本 文 将 日 尺度 CPAP 
累加 至 月 尺度 ,用 于 计算 干旱 指数 SPEI。 

1.2.2 ERAS-Land ERA5-Land 是 由 欧洲 中 期 天 气 
预报 中 心 生 产 的 第 五 代 新 气候 高 分 辩 率 再 分 析 网 
格 数据 集 ,空间 分 辩 率 为 0.1" ,和 覆盖 范围 为 90"S~ 
90°№ 2, TAX SE GE Y А 1981 年 以 来 ,表征 温 
BE .湖泊 、 积 雪 土壤 水 分 、 辐 射 和 热量 L2 НС 
流 、 风 速 气压、 降水 量 及 植被 的 50 种 变量 的 月 动态 
数据 ,具有 时 间 序 列 长 .时空 分 辩 率 高 .变量 丰富 等 
优点 。ERA5-Land 在 天 气 气候 诊断 .遥感 数据 同 
ik TF ОГИ P 3E ERIT, ERAS- 
Land 提供 了 小 时 尺度 和 月 尺度 2 种 时 间 分 辩 率 版 
本 ,本 文 从 其 官方 网 站 (https://eds.climate.coperni- 
cus.eu/) 下 载 得 到 1981 年 1 月 至 2020 年 12 月 ERA5- 
Land 的 月 尺度 莱 散 发 (PET) 网 格 数据 用 于 计算 
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图 1 研究 区 地 形 分 布 


Fig.1 Spatial distribution of terrain in study area 


SPEI 干 旱 指 数 。 该 数据 集中 的 PET 是 基于 Penman- 
Monteith 方 程 计算 ,经 验证 其 适合 于 中 国 的 干旱 监测 中 。 


2 研究 方法 


2.1 标准 化 降水 蒸 散 指数 SPEI 

不 同时 间 不同 地 区 的 降水 量 差异 巨大 ,因此 
很 难 通过 降水 量 的 大 小 进行 时 空 尺度 的 比较 。 
标准 化 降水 燕 散 指数 (SPED 同 时 考虑 了 温度 和 降 
水 , 且 具 有 多 时 间 尺 度 特征 和 标准 化 特征 ,在 实际 
应 用 中 可 消除 降水 时 空 差异 ,适用 于 干旱 时 空 变化 
规律 研究 。 前 人 研究 表明 "1 ,3 个 月 时 间 尺 度 的 
SPEI 能 够 反映 干 湿 变 化 的 季节 性 特征 , 既 不 会 因为 
时 间 太 短 导致 SPEI 值 震荡 剧烈 ,也 不 会 因为 时 间 太 
长 而 导致 SPEI 值 变化 过 于 平缓 。 本 文 基 于 1981 一 
2020 4E (40 a) 的 降水 观测 插值 数据 和 ERA5-Land 
数据 集中 的 潜在 蔡 散 发 数据 ,利用 RR 语言 中 SPEI 
包 , 计 算得 到 3 个 月 时 间 尺 度 的 SPEI( 即 SPEI3 ) , 计 
算 步 又 参考 表 1。 


22 基于 三 维 聚 类 算法 的 干旱 事件 识别 方法 

本 研究 选择 冰 值 -1 作为 干旱 状态 的 识别 标准 ， 
即 当 SPEI3<-1 时 发 生 干 旱 则 代表 该 地 区 处 于 干旱 
状态 站。 基于 Guo 等 中 的 改进 型 三 维 聚 类 干旱 事件 
识别 方法 识别 干旱 事件 并 定量 化 描述 事件 的 时 空 
演化 特征 。 该 方法 具体 步骤 如 下 : 

(D 构建 三 维 数据 空间 。 以 研究 区 的 历年 逐 月 
SPEI3 为 基本 数据 ,以 经 度 (Lon) 、 纬 度 (Lat) 和 时 间 
(为 3 个 坐标 轴 ( 如 图 2 步骤 1) ,建立 SPEI3 的 三 维 
数据 空间 。 此 时 ,该 三 维 空间 的 大 小 是 NaxNwx 
№. Ж №... Nin 和 ,分别 代 表 SPEI 三 维 数据 空间 
中 纬度 格 网 数 经 度 格 网 数 和 时 间 范 围 (月 数 )。 时 
间 范 围 N, 即 为 SPEI 干 旱 指 数 计算 的 时 间 范 围 , 本 文 
时 间 范 围 N, 范 围 设 定 为 1981 年 1 月 至 2020 年 12 月 。 

(2) 干旱 图 斑 识 别 。 对 每 个 月 的 SPEI3 数据 进 
行 干旱 状态 识别 , 当 某 格 点 SPEI3 的 值 <-1 时 则 认 
为 是 干旱 状态 并 将 其 标记 为 干旱 格 点 。 然 后 ,识别 
当前 月 所 有 干旱 格 点 并 将 相 邻 干旱 格 点 标记 为 一 
个 干旱 图 斑 。 
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表 1 SPEI 计 算 步 又 
Tab.1 The calculation steps of SPEI 


步骤 公式 备注 
1 D,- P, - PET, /表示 月 份 ; D 表示 第 /月 的 降水 ( P ) 与 该 月 潜在 燕 散 发 
(PET,) 的 差 值 
k е Hk ЗЕ E lan 个 D гу nr 
2 X È Rosy be X 3n WeK CP, SWERE PETE k SH BERN BE E 
1-81 e 第 i 年 月 的 累加 值 
aus Ур. ШЕ 
3 Ж i F(X) 是 根据 log-logistic 分 布 得 出 的 D, 序列 的 概率 分 布 函 数 ， 
а (55) 其 中 a.B 和 分 别 代表 尺度 系数 .形状 系数 和 origin 参数 
4 B= 2w, — wo шу ~ w, ~ w, 表示 log-logistie 分 布 的 3 个 矩 比 参数 ,采用 线性 


~ 6w,- w,- 6w, 


(w= 2w.) 
Ca 
netta) 
5 p=1- F(x) 


J=2ln(p) , p<0.5 
v= 
-2In( -p» P? 05 


аш 


2 
Cy T Cw t cw 


ETIE (L-moment HAZ ИН; т Z& Gamma 21i В 


Pp 是 标准 化 的 分 布 函 数值 ;w 是 一 个 参数 ; со 22.515517, сұ- 
0.802853, с, -0.010328, d, -1.432788, а, =0.189269 和 а,- 
0.001308 ,相关 参数 来 自 于 Vicente-Serrano 5? 


ІНІН (0) 
А 


t N Р 


P» 


纬度 (Lat) 


SPEI =w- 一 一 一 一 一 一 一 
1l+diwt+d,w + dsw 
步 又 1 n 
ГА H 
EFEO - 
BE 
t 
7 
же 
pA 
Z 
Р тур % 
纬度 (Lat) 
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图 2 基于 三 维 聚 类 算法 的 干旱 事件 识别 步骤 
Fig.2 Drought identification steps based on the three dimensional clutering algorithm 


(3) 干旱 图 斑 筛 选 。 在 实际 干旱 事件 分 析 过 程 
中 , 履 盖 面 广 .持续 时 间 长 的 干旱 是 研究 的 重点 , 干 
旱 碎 斑 影 响 有 限 且 会 扰乱 干旱 事件 的 识别 。 为 避 
免 王 旱 碎 斑 对 干旱 事件 识别 和 特征 定量 化 分 析 造 
成 的 影响 汪 ,本 文 仅 考 虑 面积 大 于 整个 研究 区 面积 
191.6%(15.36х10 km") H FERIRE, 

(4) 干旱 事件 识别 。 三 维 干旱 事件 是 在 干旱 指 
标 三 维 空间 中 ,聚集 在 一 起 相互 连通 的 干旱 指标 体 


素 的 集合 。 为 确定 一 个 干旱 事件 ,需要 定义 相 邻 事 
件 干旱 图 斑 的 重 肥 面积。 根据 Guo 等 的 研究 , 当 
两 个 相 邻 图 层 中 ,干旱 图 斑 的 重 县 面积 >6400 km* 
时 ,将 其 认定 为 同一 干旱 事件 。 从 第 2 个 月 开始 重 
复 步 又 2 直到 所 有 月 份 的 图 层 检查 完成 ,提取 干 早 
体 素 并 赋予 其 唯一 编号 , 即 一 场 三 维 干旱 事件 。 
23 干旱 事件 特征 变量 

本 研究 基于 以 下 6 种 干旱 特征 来 度量 干旱 
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事件 : 

(1) 干旱 持续 时 间 (Drought Duration, рр ) 是 指 
处 于 干旱 状态 下 的 时 长 , 即 干旱 开始 时 间 和 干旱 结 
束 时 间 之 间 的 月 份 数 。 

(2) FEJE. (Drought Severity,DS ) 是 指 
干旱 事件 期 间 SPEI 累加 值 的 绝对 值 , 公 式 如 下 : 


DD 
DS = > ISPEII (1) 
ial 


式 中 :DD 表 示 某 一 干旱 事件 的 持续 时 间 ;SPEL 表 示 
第 ;个 月 的 SPEI 值 。 

(3) FE ZIE (Drought Intensity, D1) 是 指 干 
旱 期 间 SPFI 值 的 平均 值 ,是 干旱 严重 度 DS 与 干旱 
持续 时 间 DD 的 比值 ,公式 如 下 : 


- ізі -DD 
DI- DD DD (2) 


式 中 :DD 表 示 某 一 干旱 事件 的 持续 时 间 ;DS 表 示 某 
一 干旱 严重 度 ;SPEL 表 示 第 ;个 月 的 SPEI 值 。 

(4) FEH (Drought Area, DA ) 又 称 干旱 范围 
(Drought Spatial Extent ) 或 者 干旱 影响 范围 (Drought 
Influenced Area) ,是 指 某 次 干旱 事件 所 涉及 的 最 大 
范围 , 即 干旱 事件 期 间 曾 处 于 干旱 状态 的 最 大 面积 。 

(5) 干旱 质心 (Drought Centroids , DC) 为 SPEI 
值 在 三 维 空间 中 的 加 权重 心 , 它 对 于 研究 干旱 的 时 
空 分 布 规律 有 着 重要 的 作用 。 为 刻画 干旱 事件 的 
发 展 迁 移 过 程 , 同 时 计算 逐 月 干旱 质心 。 

(6) 干旱 迁移 路 径 长 度 (Drought track length, 
DL) ”为 干旱 轨迹 的 长 度 , 即 逐 月 干旱 质心 移动 
距离 的 总 和 ,基于 Haversine 公式 Chttps;//en.wikipe- 
dia.org/wiki/Haversine_formula ) 计 算 两 个 质心 之 间 的 
距离 d, 其 计算 过 程 见 如 下 公式 : 
d= 


27 [A 一 入 
СЕЛІ сов Ф," a E j 


(3) 
DL- Уа (4) 


式 中 :d 表 示 地 球 上 两 点 之 间 的 距离 ;r 表 示 地 球 的 
半径 ; o, 和 gq, 分 别 表示 点 1 的 纬度 和 点 2 的 纬度 ; 
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3 结果 与 分 析 


3.1 基于 三 维 聚 类 算法 的 干旱 识别 结果 

基于 SPFI3 ,通过 三 维 聚 类 算法 的 干旱 识别 方 
法 ,在 1981 年 到 2020 年 中 国 地 区 共识 别 出 持 续 2 个 
月 及 以 上 的 干旱 事件 102 场 。 表 2 列 出 了 干旱 严重 
度 最 高 的 20 场 干旱 事件 及 其 时 空 特征 ,定量 描述 了 
每 一 场 干旱 事件 的 发 生 时 间 结束 时 间 、 严 重 程 度 、 
影响 范围 .质心 以 及 质心 路 径 长 度 。 经 统计 ,在 102 
场 干旱 事件 中 ,有 12.7% 的 干旱 历时 长 达 1a 以 上 ， 
10.8% 的 干旱 事件 覆盖 面积 超过 480x10'km2 , 占 我 
国 国土 面积 的 50%。 同 时 ,历时 超过 12 个 月 的 干旱 
事件 的 干旱 严重 度 均 位 于 所 有 干旱 事件 的 前 10% ， 
这 表明 干旱 历时 是 导致 干旱 严重 性 的 重要 因素 。 
综 上 可 知 , 近 40 a 来 中 国 的 干旱 事件 在 时 间 上 具有 
连年 性 和 在 空间 上 存在 普遍 性 。 
32 干旱 特征 变化 分 析 

图 3 显示 了 102 场 干旱 事件 质心 的 空间 分 布 ， 
图 中 圆 形 的 大 小 和 颜色 分 别 代表 干旱 事件 的 严重 
度 和 历时 。 干旱 事件 的 质心 在 研究 区 内 均匀 分 布 ， 
长 历时 、 强 烈度 的 大 规模 干旱 事件 主要 集中 在 青 
海 .甘肃 .内 蒙古 等 中 国 中 部 地 区 ,这 种 "中 间 区 域 
效应 ”主要 是 由 研究 区 域 边 界限 制造 成 的 2 ,小 规 
模 干 旱 事件 多 集中 在 南部 地 区 。 

根据 干旱 事件 中 干旱 质心 移动 的 方向 ,绘制 
1981 一 2020 年 我 国 102 场 干旱 事件 历时 -频次 的 空 
间 分 布 ( 图 4) ,其 中 不 同 的 颜色 表示 干旱 历时 。 从 
堆积 柱状 图 和 玫瑰 图 来 看 ,中国 发 生 的 干旱 更 倾向 
于 自 东 向 西 发 展 。 不 同方 向 上 的 干旱 事件 持续 时 
间 并 无 明显 规律 ,在 东 西南 . 北 、 东 北 西北 、 东 南 
和 西南 8 个 方向 中 ,西部 地 区 发 生 的 干旱 事件 累计 
值 最 高 ,达到 36 场 , 占 干 旱 总 次 数 的 35.3%, 若 加 之 
西北 西南 方向 的 干旱 事件 ,共计 有 干旱 事件 64 场 ， 
占 干 旱 总 次 数 的 62.7%。 其 中 西北 向 干旱 事件 有 16 
场 ,西南 向 干旱 事件 有 12 场 。 其 他 方向 的 干旱 事件 
略 显 弱势 , 东 、 东 北 、 北 、 南 、 东 南 的 干旱 次 数 分 别 为 
12.8、10.2 场 和 6 场 。 

根据 每 场 干旱 事件 的 逐 月 干旱 严重 度 、 逐 月 干 


A 入 分别 表示 点 1 的 经 度 和 点 2 的 经 度 ;DD 表 示 
某 一 干旱 事件 的 持续 时 间 ;d, 表 示 1 月 的 质心 与 1+1 
月 的 质心 之 间 的 距离 。 


旱 面 积 、 逐 月 干旱 烈度 ,绘制 1981 一 2020 年 干旱 时 
间 的 时 间 演 变 情况 ,通过 不 同 的 颜色 标记 不 同 历时 
的 干旱 事件 ,和 矩 形 表 示 单 次 干旱 事件 所 影响 的 面 
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图 3 干旱 质心 空间 分 布 
Fig.3 Spatial distribution of drought centroids 
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图 4 1981 一 2020 年 中 国 干旱 事件 动态 特征 干旱 历时 (a) 与 干旱 移动 方向 (b) 
Fig. 4 Drought dynamic characteristics (a) the stacked column of drought duration (b) the directional rose diagram 
of drought direction іп mainland China during 1981-2020 


耻 。 如 图 5 所 示 , 从 时 间 分 布 上 看 ,各 个 干旱 事件 的 
时 间 重 又 度 较 高 , 且 不 同 地 区 同时 发 生 干 旱 事件 的 
情况 较 多 。 从 长 时 间 序 列 的 分 布 密度 来 看 ,覆盖 范 
围 广 旦 严重 度 高 的 干旱 事件 较 多 的 发 生 于 2005 一 
2010 年 之 间 。 干 旱 严重 度 与 干旱 面积 的 相关 性 较 
高 ,但 烈度 与 历时 和 严重 度 的 相关 性 较 弱 。 这 说 明 
干旱 事件 发 生 和 覆盖 面积 越 广 更 易 引 发 高 严重 度 干 
旱 事件 。 比 如 ,2008 年 8 月 至 2009 年 9 月 发 生 的 干 


旱 事件 ,干旱 面积 覆盖 了 研究 区 80% 的 地 区 ,其 受 
旱 面 积 大 ,局 部 灾情 重 等 特点 ,使 其 成 为 研究 区 
ir 40 a 来 最 严重 的 干旱 事件 。 根 据 干旱 面积 \ 干 旱 
严重 度 及 干旱 烈度 时 间 演 化 情况 来 看 (图 5), 在 102 
场 干 旱 事件 中 , 约 40% 的 干旱 事件 呈现 出 干旱 面积 
和 干旱 严重 度 的 多 峰 现 象 , 即 一 场 王 旱 事件 过 程 
中 ,其 干旱 面积 和 严重 度 变化 出 现 多 个 峰值 。 超 过 
半数 (57.84% ) 的 干旱 事件 在 干旱 烈度 变化 中 呈现 
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图 5 1981 一 2020 年 中 国 和 干旱 事件 的 严重 度 .面积 .烈度 的 时 间 演 变 


Fig. 5 Temporal evolution of the severity, area апа intensity of drought events in China, 1981-2020 


多 峰 现 象 。 不 同 特征 的 干旱 事件 多 峰 现 象 不 利于 
二 时 事件 发 展 的 预 佑 ,这 是 由 于 当 某 一 干旱 特征 
(干旱 烈度 、 干 旱 严 重度 、 干 旱 面积 ) 由 升 转 降 时 并 
不 能 代表 干旱 的 结束 ,也 许 会 产生 新 的 上 升 趋势 。 
受气 候 变 化 的 影响 ,中 国 地 区 三 维 干旱 特征 变 


形 为 单个 干旱 事件 的 影响 面积 。 


化 情况 较为 复杂 ,图 6 以 5 a 为 周期 统计 了 干旱 事件 
的 累计 频次 .历时 均值 .烈度 均值 .严重 度 均值 和 面 
积 均值 。 可 以 看 出 ,40 a 来 , 单 场 干旱 事件 的 持续 月 
数 均值 随时 间 变 化 呈现 明显 的 抬升 趋势 ;干旱 烈度 
均值 稳定 在 0.6 左 右 ; 干 旱 历 时 均值 .干旱 严重 度 均 
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图 6 1981—2020 4 rp p 
Fig. 6 Decadal change plot of drought characteristics in China during 1981-2020 


f& .干旱 面积 均值 趋势 大 致 相同 , 呈 先 抬升 后 轻微 


司 干旱 特征 年 代 际 变化 


10 月 以 后 逐渐 加 重 ,2009 年 1 月 下 旬 各 地 旱情 达到 


下 降 的 形势 。 综 合 来 看 ,2000 年 之 后 长 历时 干旱 事 
件 发 生 次 数 明显 增多 ,2005 一 2010 年 间 发 生长 历时 
干旱 事件 16 场 ,平均 干旱 历时 和 面积 分 别 为 6.7 个 
月 和 2.8x10' kn, 历时 长 .烈度 高 .影响 面积 广 的 干 
旱 事件 大 多 发 生 在 2005 一 2010 年 。 
3.3 典型 三 维 干旱 事件 分 析 

为 分 析 典 型 干旱 事件 的 空间 格局 和 时 间 滨 变 
特征 ,本 文 根 据 三 维 干旱 特征 值 的 不 同 , 从 最 严重 
的 20 场 干旱 事件 中 ,再 次 选取 出 严重 度 最 高 .烈度 
最 高 .历时 最 长 \ 最 短 .同等 严重 度 下 历时 较 短 .同等 
严重 度 下 烈度 较 强 的 6 场 典 型 干旱 事件 ,并 在 表 2 
中 以 淡 黄色 背景 标记 。 图 7 和 图 8 展示 了 典型 干旱 
事件 三 维 结构 及 特征 变量 的 时 间 趋 势 。 由 图 8 可 
知 ,在 单 场 干 旱 事件 中 ,干旱 严重 度 .干旱 面积 和 干 
旱 烈 度 的 逐 月 变化 趋势 基本 一 致 。 

通过 干旱 事件 d( 图 7d) 可 以 明显 看 出 ,作为 研 
究 时 段 内 最 严重 的 干旱 事件 ,本 场 干旱 事件 的 干旱 
强度 高 .影响 面积 广 ,从 中 国 东北 地 区 逐步 发 展 至 
中 国 西部 地 区 ,而 后 其 干旱 质心 向 东 移 动 ,最终 在 


了 最 严重 的 程度 。 

干旱 事件 {( 图 790) 是 近 40 a 来 发 生 的 历时 最 长 
的 干旱 事件 , 共 历 时 22 个 月 。 在 历时 较 长 的 干旱 事 
件 中 ,干旱 严重 度 .干旱 面积 等 特征 往往 会 产生 多 
个 峰值 。 这 种 多 峰 现 象 不 利于 干旱 发 生发 展 过 程 
的 预测 ,比如 干旱 严重 度 等 特征 的 暂时 降低 可 能 被 
错误 地 认为 是 干旱 进入 消亡 期 。 此 次 干旱 事件 主 
要 发 生 在 中 东部 地 区 ,累计 覆盖 面积 为 3.78x105 
km , 远 小 于 干旱 事件 d 的 面积 。 在 中 国 应 急 管 理 部 
发 布 自然 灾害 基本 情况 中 描述 到 ,2019 年 东北 地 区 
遭遇 春 旱 ,山西 .河南 等 地 出 现 阶段 性 夏 伏旱 , 华 
北 、 西 南 和 东北 地 区 夏 伏旱 人 3, 其 描述 的 干旱 影响 
范围 与 本 文 识别 的 干旱 面积 基本 一 致 。 

于 旱 事件 a( 图 7a) 的 历时 与 干旱 事件 4( 网 7d) 
相同 ,但 其 干旱 严重 程度 远 小 于 干旱 事件 4( 图 7f)， 
究 其 原因 ,主要 是 其 干旱 面积 不 同 。 干 旱 事 件 a 从 
中 国 中 部 地 区 逐步 向 西 发 展 ,干旱 受灾 范围 相对 集 
中 ,其 峰值 出 现在 本 场 干旱 发 生 的 第 8 个 月 ,其 单 月 
干旱 影响 面积 达到 了 1.54x105km2。 由 此 可 知 , 在 干 


中 国 中 部 地 区 消失 ,覆盖 面积 达 8.00x10'knm2 , 约 占 
中 国 陆 地 面积 的 80% ,历时 长 达 14 个 月 。 在 2009 
年 2 月 ,此 次 干旱 事件 达到 覆盖 面积 上 的 峰值 ,其 覆 
盖 面 积 超过 5.00x10'km*。 在 前 人 人 研究 中 验证 了 
2008 年 10 月 至 2009 年 3 月 我 国 北方 地 区 出 现 的 秋 
冬 连 旱 现 象 汪 ” ,人 研究 中 描述 了 此 次 旱情 在 2008 年 


旱 事件 历时 相同 的 情况 下 ,干旱 禾 盖 面积 都 对 干旱 
事件 的 严重 程度 起 决定 性 作用 。 干 旱 事 件 c( 图 
7e) ,发 生 在 中 国 的 东北 地 区 ,与 干旱 事件 a 的 干旱 
严重 度 相同 ( 均 为 1.74x10’ km?) ,值得 注意 的 是 , 尽 
管 干旱 事件 c 的 历时 短 且 获 盖 面积 小 ,但 其 干旱 烈 
度 却 比 干 旱 事件 a 高 出 了 0.58, 说 明 在 干旱 严重 度 
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相同 的 状态 下 ,干旱 烈度 越 小 ,干旱 面积 则 越 大 , 反 

干旱 事件 b( 图 7b) 的 严重 度 在 所 识别 干旱 事件 
中 位 列 第 4, 与 同等 严重 度 的 干旱 事件 相 比 ,此 次 干 
早 事件 的 历时 较 短 ,表明 此 次 干旱 事件 发 生 短促 而 


烈度 最 强 干 旱 事件 (2007 年 6 月 至 2008 年 5 月 ) 的 时 
空 动态 演变 过 程 。 图 9 绘制 了 干旱 事件 发 生 至 消亡 
的 时 空 动 态 演变 全 过 程 及 其 逐 月 迁移 路 径 , 点 代表 
每 个 月 干旱 集群 的 质心 ,连接 质心 之 间 的 箭头 为 干 
旱 迁 移 路 径 。 由 图 可 知 ,本 场 干 旱 事件 的 源头 在 内 


剧烈 。 该 干旱 事件 发 生 自 陕西 山西 以 及 内 蒙古 交 


蒙古 北部 地 区 ,发 生 干 旱 时 质心 在 赤峰 市 ,随后 干 


汇 处 ,而 后 迅速 草 延 至 全 国 ,直至 1999 年 4 月 ,干旱 


旱 迅 速 昌 延至 东北 地 区 全 域 , 干 旱 质心 向 东北 方向 


质心 移动 至 西藏 地 区 ,而 后 逐渐 消亡 。 干 旱 事件 e 
(图 7e) ,此 次 干旱 事件 的 历时 为 6 个 月 ,是 中 
1980 一 2020 年 最 严重 的 20 场 干旱 事件 中 历时 最 短 


迁移 至 白城 市 ;8 月 至 9 月 ,干旱 继续 向 北 发 展 , 干 旱 
严重 程度 持续 增加 ,旱情 最 为 严重 , 单 月 干旱 面积 
达到 峰值 (1.38x10'km”); 10 月 干旱 分 散 为 两 个 部 分 ， 


的 干旱 事件 ,但 其 干旱 烈度 处 于 第 2 位 。 与 干旱 事 
TF b 的 干旱 特征 一 致 ,此 干旱 事件 的 发 生发 展 同样 
的 短促 而 剧烈 ,值得 注意 的 是 ,干旱 事件 e 在 一 个 月 
内 消亡 ,消亡 期 十 分 短暂 。 

以 上 干旱 事件 与 有 关 人 研究 文献 “相符 合 , 每 场 
干旱 事件 的 发 生发 展 时 间 和 影响 范围 与 记载 的 一 
致 性 较 高 ,进一步 验证 了 本 文 提 出 的 三 维 聚 类 算法 
干旱 识别 是 合理 可 靠 的 。 

为 分 析 典 型 干旱 事件 的 发 生发 展 过程 逐 月 空 
间 影响 范围 等 特征 ,本 文 进一步 分 析 了 人 研究 时 段 内 


旱情 由 此 开始 衰减 ;11 月 至 次 年 1 月 ,干旱 中 心 向 西 
迁移 至 齐齐哈尔 市 ;2 月 至 3 月 ,干旱 中 心 向 东 渐 渐 
移动 ;4 月 至 5 月 ,干旱 中 心 持 续 向 北 移动 ,最 终 移 动 
至 呼伦贝尔 市 北部 ,并 最 终 在 此 消亡 。 综 合 来 看 ， 
本 次 干旱 事件 从 2007 年 6 月 至 次 年 5 月 共 历 时 12 
个 月 ,此 干旱 生长 期 较 短 (2007 年 6 月 至 2007 年 9 
月 , 共 4 个 月 ) ,衰减 时 间 较 长 (2007 年 10 月 至 2008 
年 5 月 , 共 8 个 月 ), 旱 情 经 历 了 生长 期 -强化 期 - 峰 
值 期 -衰减 期 -消亡 期 5 个 阶段 ;干旱 质心 移动 方向 
由 南 向 北 迁 移 ,迁移 路 径 大 致 为 赤峰 市 一 白城 市 一 
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图 7 典型 干旱 事件 三 维 结构 


Fig.7 The three dimensional structure of typical drought events 


202308.00251v1 


ChinaXiv 


月 干旱 面积 /105km? 月 于 旱 面积 /105km? 
月 干旱 严重 度 /105 km? 月 干旱 严重 度 /105 km? 


月 干旱 面积 /105km? 
月 干旱 严重 度 /105 km? 


(а) 2007-06 


(е) 2007-10 


3 


№ 


N 
RUD 
月 干旱 烈度 
H ЛО kar? 


о N ы о о 


12H 1H 2H 3 月 


6H 8H 


昌 潇 雨 等 :基于 三 维 识别 的 中 国 干 旱 事件 演变 特征 分 析 


ЕШ 月 干旱 面积 
(а) 1986-08--1987-09 


一 一 月 干旱 严重 度 —— 月 干旱 烈度 


8H 10H 12H 2H 4H 6H 8H 


(c) 2007-06--2008-05 
T T 


10H 12H 
(е) 2012-12-22013-05 


2H 4H 


月 干旱 烈度 
月 干旱 面积 /105 km? 


图 8 典型 干旱 9 


月 干旱 烈度 
月 干旱 面积 /105 km? 


4H 5H 


(b) 1998-10--1999-06 


ChinaXiv 合 作 期 刊 


859 


ші 12 


E 6 MM 18 
5 4 у> 1.6 E: 
p 
H 2 14 
E 0 12 

108118 12H1H 2 月 3H 4H 5H 6H 
a (4) 2008-08--2009-09 
Б 10 T——1——1 T——1 T 1.7 
FEE Е 
EN wc 
І Шш 
T o 
um 

8H 10H 12H 2H 4H 6H 8 月 

" (f) 2018-10— 2020-07 
E 4 2.0 
Š 3 18 $ 
Шо; 16 € 
EN H 
ПЕ 1 1.4 um 
H 
um 


10H 1H 4H ТА 10H ІЛ 4H 7H 


和 件 特 征 变 量 的 时 间 变 化 


Fig. 8 Temporal changes of the characteristics of typical drought events 
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和 件 时 空 动态 演变 过 程 及 其 逐 月 迁移 路 径 


Fig.9 Evolution of the spatial and temporal dynamics of drought events from June 2007 to May 2008 and 


尔 市 北部 。 此 次 干旱 事件 与 中 国 气 象 干旱 图 集 “ 


their migration paths from month to month 


大 庆 市 一 呼伦贝尔 市 南部 一 齐齐哈尔 市 一 呼 伦 由 


中 2007 年 秋季 干旱 、2008 年 春季 干旱 重合 度 较 高 ， 
进一步 验证 了 三 维 聚 类 算法 识别 干旱 事件 的 可 


THE 


4 


结论 与 讨论 


本 文 以 中 国 为 研究 区 ,从 三 维 视角 出 发 ,识别 
1981 一 2020 年 间 的 干旱 事件 ,并 系统 分 析 干 旱 事件 
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的 历时 .严重 度 .烈度 .影响 面积 等 基本 特征 和 移动 
路 径 .发 展 方向 等 动态 特征 。 基 于 三 维 聚 类 算法 的 
干旱 事件 识别 方法 能 够 为 揭示 干旱 时 空 变化 规律 
提供 新 途径 ,研究 结论 有 助 于 深化 中 国 干旱 事件 发 
生发 展 过 程 的 认识 和 理解 。 本 文 主要 结论 如 下 : 

与 历史 文献 资料 对 比分 析 ,结果 表明 三 维 聚 类 
算法 能 有 效 识 别 区 域 气象 干旱 事件 ,1981 一 2020 年 
间 共 识别 历时 超过 2 个 月 的 干旱 事件 102 场 。 研 究 
时 段 内 ,最 严重 干旱 事件 发 生 于 2008 年 8 月 至 2009 
年 9 月 期 间 , 共 历 时 14 个 月 ,影响 了 中 国 超 80% 的 
面积 。 不 同 干旱 事件 的 时 空 特征 差异 显著 ,同一 个 
干旱 事件 逐 月 烈度 、 逐 月 严重 度 等 特征 随时 间 的 变 
化 区 别 也 较 明 显 。 

基于 干旱 事件 的 质心 分 布 情况 以 及 历时 图 、 移 
动 方 向 玫瑰 图 分 析 发 现 , 中 国 发 生 的 长 历时 、 强 烈 
度 的 大 规模 干旱 事件 主要 集中 于 中 西部 地 区 , 且 
62.7% 的 干旱 事件 更 倾向 于 自 东 向 西 发 展 。 通 过 逐 
月 分 析 干 旱 事件 的 特征 变化 ,发 现 中 国 发 生 的 干旱 
事件 时 间 重 又 度 较 高 。 从 长 时 间 序 列 的 分 布 密 度 
来 看 , 窗 盖 范围 广 量 严 重度 高 的 干旱 事件 较 多 的 发 
生 于 2005 一 2010 年 之 间 。 在 102 场 干旱 事件 中 , 约 
40% 的 干旱 事件 呈现 出 干旱 面积 和 干旱 严重 度 的 多 
EMA ,超过 半数 (57.84% ) 的 干旱 事件 在 干旱 烈度 
变化 中 呈现 多 峰 现 象 ,干旱 事件 的 多 峰 现象 不 利于 
于 旱 事 件 发 展 的 预 估 。 

本 文 针 对 干旱 过 程 在 时 空 维度 下 的 演变 特征 
开展 研究 仍 存 在 一 些 待 改 进 之 处 。 首 先 ,SPEI 干 早 
指数 是 基于 降水 和 潜在 蒸 散发 计算 得 到 ,本 质 上 反 
IU Zi РА ,并 未 考虑 下 垫 面 植被 类 型 .地 
形 等 因素 的 影响 。 其 次 ,部 分 研究 ”表明 ,以 研究 
区 1.6% 的 面积 比重 作为 最 小 干旱 聚 类 面积 阔 值 能 
够 有 效 提 取 三 维 干旱 事件 并 避免 干旱 聚 类 碎 斑 出 
现 , 然 而 不 同 立 值 的 大 小 会 导致 干旱 事件 提取 结 
的 差异 ,其 阔 值 选取 的 合理 性 仍 有 等 进一步 验证 。 
最 后 ,三 维 聚 类 干旱 识别 算法 虽然 能 够 有 效 避 人 免 干 
旱 事 件 的 动态 信息 丢失 的 问题 ,但 同时 面临 算法 
相对 较 复杂 的 问题 。 
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Characterization of the evolution of drought events in China 
based on 3D identification 


LYU Xiaoyu, GUO Hao, MENG Xiangchen, ВАО Anming, TIAN Yunfei, ZHU Li 
(1. School of Geography and Tourism, Qufu Normal University, Rizhao 276826, Shandong, China; 
2. Xinjing Institute of Ecology and Geography, CAS, Urumqi 830011, Xinjiang, China) 


Abstract: Globally, China is one of the country' s most frequently and severely affected by drought. These 
frequent drought events have subsequently caused serious economic, social development, and ecological 
environment losses. Drought simultaneously leads to alterations in both space and time, and this paper aims to 
identify drought events and quantify their spatial and temporal dynamic evolution using a three- dimensional 
clustering algorithm for mainland China, from 1981-2020. The three- dimensional clustering algorithm can be 
used to effectively identify drought events and describe their dynamic processes. From 1981-2020 there were 
102 drought events lasting 2 months or more in mainland China. Spatially, the trajectory of drought events was 
found to show a tendency to develop from east to west. Temporally, there are high time overlaps between 
different drought events. In addition, drought events with wide coverage and high severity were concentrated in 
the period from 2005-2010. The findings of this paper will help to elucidate the spatial and temporal evolution 
patterns of drought events and provide scientific references for drought monitoring and drought risk management 
in mainland China. 

Keywords: three-dimensional clustering algorithm; drought event; three-dimensional evolutionary features; dy- 
namic evolution; SPEI; China 


